From Evernote: |
TJGR Using gill expression get interval file for MBD Gill |
1) Run tophat in iPlant {Gill Expression Zhang et al. v9_90}
SELECT Column1,Column4,Column5,Column14 FROM [sr320@washington.edu].[table_TJGR_Gil_Exp_TH_coverage_CDS.csv]
SELECT * FROM [sr320@washington.edu].[Gil_Exp_coverage_CDS_SummaryStats]
SELECT Column9, sum(Column10) FROM [sr320@washington.edu].[table_TJGR_CDS_intersect_gillMBDmCpG]
SELECT * FROM [sr320@washington.edu].[Gill_Done]
Select "Gene_name" From [Gill_Methylome]
Select * From [sr320@washington.edu].[Gill_Done_2]
Select * From [sr320@washington.edu].[Gill_Done_2]
--
converting oyster CDS gff to gtf for tophat?
Started adding GTF file that has CDS and mRNA data - looks for splice junctions.
Results
------
2) Get coverage on CDS
-split
Default Output:
After each entry in B, reports:
1) The number of features in A that overlapped the B interval.
2) The number of bases in B that had non-zero coverage.
3) The length of the entry in B.
4) The fraction of bases in B that had non-zero coverage.
cleaned up in excel
column on end - num A overlap / length of B
Where Column14 > 2
Now using Galaxy to filter "high expression"
Using last column(column14) - 0.7 is about 22% CDS with highest expression
format: gff, database: oyster_v9_90
Filtering with c14>=0.7, kept 22.65% of 196691 valid lines (196691 total lines).
Back in SQL
SELECT Column9, sum(Column14),count(Column14),(stdev(column14)/avg(Column14)),var(column14),avg(Column14),stdev(column14),min(column14),max(Column14) FROM [sr320@washington.edu].[table_TJGR_Gil_Exp_TH_coverage_CDS.csv]
Group by Column9
Where Column4 ]] > 0
avg Expression > 0
NOW have Gill Expression coverage on CDS with Stats; including CV
Once have mCpG data for gill MBD library, will see if there is a relationship
./intersectBed -c -a /Users/sr320/Desktop/oyster.v9.glean.final.rename.mRNA.gff -b /Users/sr320/Desktop/TJGR_GillMBD_mCpG_10xcov.bed > /Users/sr320/Desktop/TJGR_mRNA_intersect_gillMBDmCpG
Group by Column9
Join mCpG data will all CG - mRNA
JOINING
[sr320@washington.edu].[TJGR_Gill_meth_CDSexpression_gene]
--
SELECT GENEID,Gil,swissprot
FROM [sr320@washington.edu].[table_Table S14.csv]
Join with Gill Expression
CDS count = 10
Where Percent_Meth > 20
and CG > 10
and express < 10
Where Percent_Meth < 40
and CG ]] > 10
and CDScount = 10
Where Percent_Meth ]] > 60
and CG ]] > 10
and CDScount = 10
---
DH BB3 Expression
---
3) Get extraCDS windows in mRNA (maybe flank)
3b) Get coverage on extraCDS windows
----
Join in SQLshare and select datasets with similar overall expression, x # exons,
stats to try to identify alternative splicing…
--
What about relationship with SNPs
#get that.
----------
From: Steven Roberts <sr320@u.washington.edu>
Date: Mon, Feb 18, 2013 at 11:54 AM
To: Mackenzie Gavery <mgavery@uw.edu>
Steven Roberts